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Abstract. In a companion paper, we established nonlinear stability with detailed diffusive rates 
of decay of spectrally stable periodic traveling-wave solutions of reaction diffusion systems under 
small perturbations consisting of a nonlocalized modulation plus a localized perturbation. Here, 
we determine time-asymptotic behavior under such perturbations, showing that solutions consist 
to leading order of a modulation whose parameter evolution is governed by an associated Whitham 
averaged equation. 



1. Introduction 

In this, the second part of a two-part series of papers, we determine the time-asymptotic behav- 
ior of spectrally stable periodic traveling-wave solutions of reaction diffusion systems under small 
perturbations consisting of a nonlocalized modulation plus a localized perturbation, showing that 
solutions consist of an (in general nonlocalized) modulation governed by the formal second-order 
Whitham averaged equations plus a faster-decaying localized residual. 

In the first part, [JNRZl], we established nonlinear stability under such perturbations, together 
with detailed rates of decay, using a refined version of the argument used in [JZl] to show stability 
under localized perturbations. In each of these cases, the basic approach was to separate out from 
the linearized solution operator slowly-decaying terms corresponding to linear modulations, then 
estimate separately the modulational and nonmodulational parts of the solution. 

Here, we show that, using the same basic linear estimates as in [JNRZl], but a one-order- 
higher decomposition of the solution operator, we may obtain also asymptotic behavior in terms 
of a modulation determined implicitly in terms of a system of integral equations. By a further 
linear estimate in the spirit of [JNRZl], we find that the principal part of the local wavenumber 
appearing in this modulation obeys a convected Burgers equation that is asymptotically equivalent 
to the Whitham modulation equations (in the sense of [SI, DSSS, W]) . 

This not only gives rigorous validation of the formal Whitham approximation in the strong sense 
of showing that it describes time-asymptotic behavior, but, through the explicit prescription of 
initial data coming from our analysis, also gives new information of predictive value not available 
from the formal asymptotic derivation. Indeed, as discussed in [BJNRZ2], the connection between 
initial perturbation and initial values for the Whitham equations has remained for a long time 
somewhat mysterious. As a side-consequence, we show that the decay rates of [JNRZl] are sharp. 

Consider a periodic traveling-wave solution u{x, t) = u{ki,{x — ct)) of reaction diffusion equation 
ut = Uxx + f{'^)^ or, equivalently, a standing-wave solution u{x,t) = u{x) of 

(1.1) k^ut = kluxx + f{u) - ioux, 
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where a;(A;*) := —k^c{kt,) is the temporal frequency, c is the speed of the original traveling wave, 
and wave number A;* is chosen so that 

(1.2) u{x + l) = u{x). 

Here and throughout the paper, all periodic functions are assumed to be periodic of period one. 
We make the following standard genericity assumptions: 
(HI) / e C^(R) for some K > 3. 

(H2) Up to translation, the set of 1-periodic solutions of (1.1) (with k replacing A;*) in the vicinity 
of u, k = A;^, forms a smooth 1-dimensional manifold {u{k, •)} = {'U*^(-)}, cu = u){k). 

Linearizing (1.1) about u yields the periodic coefficient equation 

(1.3) k^vt = ki^Lv := {kldl - Lo{ki,)dx + b)v, b{x) := df{u{x)). 
Introducing the family of Floquet operators 

(1.4) k,L^ := e-*«"A,Le^«" = kl{d., + i^)' " ^{hW, + iO + b, 
operating on periodic functions on [0, 1], determined by the defining relation 

(1.5) Lie'^""/) = e^«^(Lg/) for / periodic, 

we define following [SI, S2] the standard diffusive spectral stability conditions: 
(Dl) = is a simple eigenvalue of Lq. 
(D2) a{L^) C {A I 5RA < for some constant 9 > 0. 

The following stability result was established in [JNRZl], generalizing results of [SI, S2, JZl]. 
Proposition 1.1 ([JNRZl]). Let K > 3. Assuming (H1)-(H2) and (D1)-(D2), let 

Eq •— \\uq{- - /lo(-)) - '"(■)llLi(M)ni?^(M) + ll^a:^o||Li(]R)nif^(IR) 

be sufficiently small, for some choice of modulation ho. Then, there exists a global solution u{x,t) 
of (1.1) with initial data uq and a phase function ip{x,t) such that, for t > and 2 < p < oo, 

M- - i^{;t),t) - t2(.)|U.(M) < Eo{l + t)-^(i-Vp) 

\\V.,M;t)\\LPiR)<Eo{l+tr"^^'-'fP\ 

and 

(1.7) \\u{-,t) -u{-)\\Loog^), \\llj{-,t)\\Loo(TR)<Eo. 

In particular, u is nonlinearly (boundedly) stable in L°° with respect to initial perturbations vq = 
uo-u for which \\vo\\e := ^^id^hoeL^K)nHK{M_) Eq is sufficiently small. 

Recall now the formal, Whitham equation, as derived in various contexts and to varying degrees of 
accuracy in [W, HK, Se, DSSS, NRl, NR2]. By translation-invariance of the underlying equations, 

(1.8) Lou' = 0, 

so that by (Dl) the zero-eigenspace of Lq is exactly Rangeju'}. Denote by u"'^ the left, or adjoint, 
zero eigenfunction of Lq. Fixing u, A*., introduce the parametrization 

(1.9) u''{x- P) = u{k,x- P) 
of nearby periodic traveling waves, i.e., 1-periodic solutions of 

(1.10) ujik)4 - k-'ul, - fiu'^) = 0, 
uj{k) = —kc{k), for definiteness chosen in such a way that 

(1.11) {u''''{k),dku{k))LH[o,i]) = 0- 
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(This may achieved by an appropriate translation, by the fact that (u""^, u') 7^ 0.) Then, the formal 
approximate solution of ut — Uxx — f{u) = obtained by a nonlinear WKB expansion is 

(1.12) u{x,t)^u^^''^''\^{x,t)), 

where the wave number k := satisfies the Whitham equation (viscous scalar conservation law) 

(1.13) Kt - = {d{K)Kcc)x 

or equivalently, the phase ^ satisfies its integral (viscous Hamilton-Jacobi equation) 

(1-14) *t - oj{^x) = di^x)'fxx 

where, taking (1.11) into account, 

(1.15) dik) = l + 2A;(u«<^(A;),afe^x'(fc))i2([o,i]) . 

See [NRl, NR2] for a detailed derivation of this kind of nonlinear Whitham's equation in the context 
of the Saint- Venant and Kuramoto-Sivashinsky equations. 

In our context, we look for the evolution of a localized perturbation k = k^^hx in our co-moving 
frame, and so it is enough to retain the quadratic order approximants 

(1.16) k^kt + k^q{k)x = kld{ki,)kxx, 
and 

(1.17) k^ht + q{k^hx) = kld{k^)hxx 
with 

(1.18) qik) = -iuj'ik,) + c{k,))k - ]^J'{k,)k'' . 

For, as is well known, (1.13) and (1.16) are "asymptotically equivalent" in the sense that the 
difference between solutions of (1.13) and (1.16) decays faster in all U", 1 < p < 00, than does the 
solution itself, which, for an initial perturbation with nonzero integral, decays at the rate of a heat 
kernel. See [DSSS] for a direct derivation of the quadratic Whitham's equation. 

Lemma 1.2 ( [Ka, LZ] ). Let ij > be arbitrary. Let k be a solution of (1.13) with initial datum kq 
and k be a solution of (1-16) with initial datum ko = ko{ ■ /fc*) — k^,, Eq := ||feo||Lini?3(K) sufficiently 
small. Then, setting k{x, t) = ki, + k{k-),{x — c{ki,)t) , k^^t) , 

- ^)(t)||iP(M) < E^il + 1 < p < 00; 

moreover, if mo := j'^ko and i?i := £'0 + || | • | A;o||ii(iR) is sufficiently small, then 

Mt) - k^LPiR), Mt) - k4Lp{R) > \mo\ (1 + t)-^(i-Vp) . 
Proof. See Appendix A. □ 
Our main result is as follows. 

Theorem 1.3. Let 77 > and K >3. Under the assumptions of Proposition 1.1, let k and h satisfy 

the quadratic approximants (1.16) and (1.17) of the the second-order Whitham modulation equations 
(1.13) and (1.14) with initial data k\t=o = k^^dxho, h\t=o = ho. Then, for t > 0, 2 < p < 00, and ^ 
as in Proposition 1.1, 

||n(. - ^(., t),t)- u''*(^+M;t)) ( . ) < £0 ln(2 + t) (1 + t)-'^ , 

(1-19) \\kAm - miiPiR) < Eoii + t)-i(i-vp)-i+'?, 

-h{t)\\LPiR)<Eo{l + tr'^^'-'/P^+^. 
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Remark 1.4. From estimates (1.19) (i) on u and (1.6) (ii) and (1.7) on ipx and ip, we obtain 
\\ui;t) - u'=*(i+'^^(-'*»(*(-,t))||L.(M) < ^oln(2 + t) (1 + i)"^, 2 < p < oo, 

where ^'(•, t) is the inverse oi y X{y, t) := y — il^{y, t). We insure the existence of such a map by 
keeping, for any t, ||V'(t)||L°°(M) bounded and ||'0x(^)||l°c(r) small. Since ^(x,t) — {x + il){x,t)) = 
ijj{^{x,t),t) — ip{^{x,t) — ip{'^{x,t),t),t), one may translate this into a bound 

(1.20) \\u{; t) - u^'^^^i-'i) (*(., t)) IUp(k) < Eo ln(2 + t) (1 + t)-f , 2 < p < oo, 

of the form (1.12), or degrade it into^ ||u(-,t) - ii''*(i+'^-(-'*))(- +'0(-,t))||iP(K) < Eo{l + t)-^^^-^'P\ 
obtaining thereby, for 77 > arbitrary, and 2 < p < 00, the exact Whitham comparison 

(1.21) \\u{-,t) - #*+M-,0(. + /i(.,0)|U.(ffi) < Eo{l + t)-"^^^-yp)+^. 

From Theorem 1.3 and Lemma 1.2, we see immediately that the decay rates of Proposition 1.1 
are sharp. At the same time, we give rigorous validation through (1.21) of the Whitham equation 
(1.13) in the long-time limit. However, the sharpened bound (1.20) shows that we have in fact done 
considerably more. For a corresponding result in the conservative case, see [JNRZ2]. See [DSSS] 
for validation of the Whitham equation in the sense of building a family of solutions existing on 
asymptotically-large but bounded intervals and close to a given asymptotic expansion involving a 
given solution of the Whitham solution. 

Remark 1.5. Without change, we may treat reaction-diffusion systems with a diagonal diffusion. 

More generally, all proofs go through whenever the linearized operator generates a semigroup; 
in particular, our results apply to the (sectorial) Swift-Hohenberg equations treated for localized 
perturbations in [82] . It would be interesting to extend to the multi-dimensional case as in [U] . 

Note: We have been informed by Bjorn Sandstede that similar results have been obtained by 
different means by him and collaborators [SSSU] using a nonlinear decomposition of phase and 
amplitude variables as in [DSSS]. 

2. Preliminaries 
Recall the Bloch solution formula for periodic-coefficient operators, 

(2.1) {Sm{^) ■■= (e'^^aKx) = r e^«^(e*^«5(er))(x)de,' 

J—TT 

where is as in (1.4), 

(2.2) K^,x):= J]^(C + 2J7r)e^2-^■^ 

periodic, denotes the Bloch transform of g, g{^) := ^ J^e~^^^g{x)dx the Fourier transform, and 

(2.3) g{x)= r e'^-g{C,x)d^ 



the inverse Bloch transform, or Bloch representation of g € L^. The generalized Hausdorff- Young 
inequality ||tf||LP(R) ^ II^IIl'j(5,lj'([o,i])) for g < 2 < p and | + | = 1 [JZl], together with (2.3), yields 
for any 1-periodic function /(^, •) 

(2.4) II r e*«7(^, ■)d^\\LP(m < II/IIl.(€,lp([o,i])) for g < 2 < p and - + - = 1. 

J —n P 1 

We use here, ||'0(i)||LP(H) < Eo{l + i)-2(i-i/p)+2 , for 2 < p < 00, a straightforward extension of (1.7). 
^In other words, {e*^g){^,x) = {e^^ig{^,-)){x) a consequence of (1.5). 
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By (Dl), the zero eigenfunction u' of Lq is simple, whence by standard perturbation results 
bifurcates to an eigenfunction ^(^, •), with associated left eigenfunction (f){^, •) and eigenvalue 

(2.5) x{0 = aiC-de + O{\Cf), 

where a and d are real and d > by assumption (D2) and complex symmetry A(^) = A(— each 
of (f), (p, A analytic in ^ and defined for |^| sufficiently small. 

3. Linear estimates 
Refining slightly the decomposition of [JZl, JNRZl], decompose the solution operator as 

(3.1) S{t) = RP{t) + R{t), RP{t) = {v! + k^dkud^)sP{t), 
with 

(3.2) (s^mix) = r e^^^aiOe'^i^'m, ■))m[o,i])d^, 

J—n 

and 

(Rmi^) ■■= r e^«-(l-a(C))(e^«*5(e))(^)d^ + T e'^^a{0{e'^^'mmm^)d^ 

J—n J —TT 

+ re'«-a(Oe"(«)*(</.(C,x)-0(O,a;)-C5^0(O,x))(^(O,5(O)L^([o,i])^^C, 

J— IT 

where a is a smooth cutoff supported on |^| sufficiently small, 

defined for sufficiently small ^ denotes the eigenprojection onto the eigenspace Range{^(^)} bifur- 
cating from Rangej-u'} at ^ = 0, ^ the associated left eigenfunction, and n := Id — 11*'. 

Proposition 3.1 ([JZl, JNRZl]). Under assumptions (H1)-(H2) and (D1)-(D2), for all t > 0, 
2<p<oo, 

•(l+t)-^(^-V-)-^||5lUi(M), 

|5||l2(k), 

for < n < K + 1, and for some rj > and < I + 2m, n < K + 1, 



(3.3) 



(3.5) dldrs^imi^g 



L^ffi) ~ 1 (1 + i)-§(i/2-i/p)-^^ I 



(3.6) didrR{t)d^g 



<^i^f(l + *) ^^^^ ^lbllLi(K)n/f'+2-+i(R), 

L.(E) ~ I e-^^d^g\\Hi+2^+im + (1 + t)-'^^'^'-'/P^~'\\g\\LH 



Proof. The estimates on s^ were proved in [JNRZl]. Using (2.4) together with spectral projec- 
tion/direct computation for low frequencies and standard semigroup estimates for high frequencies, 

estimates on R arc proved exactly as were the estimates on S in [JNRZl], with the observation 
that the substitution of factor (0(^) — (/>(0) — S,d^(p{0)) in R for factor {(/){£,) — 0(0)) in S introduces 
an additional factor of |^| in the estimates, hence an additional (1 + t)^^/^ factor of decay. □ 

Proposition 3.2 ([JNRZl]). Under assumptions (H1)-(H2) and (D1)-(D2), for all t > and 

2<p<oo, 

(3.7) \\didrsP{t){hou')hnR) < (1 + t)-^(^-^/^)+^-'^||5./^o||Li(R), 
when I + m > 1 or else I = m = and p = oo, and, for < I + 2m < K + 1, 

(3.8) \\dldrRmhou')\\LPiR) < {l + t)-'^^'-'^P^-'\\dM\LHR)nHi+^ru+i(j^), 
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and when t < 1 



.oQ^ \\di,driRnt)-ld){hou')\\LPiR) < Wd^hoW Li(K)ni?'+2'"+i(R)) 

\\dld^{sP{t){hou') - ho)\\Lp(j^) < ||5x^o||Li{M)nL2(iR)- 

Proof. Again, estimates on were proved in [JNRZl] whereas the proof of those on and R 
goes exactly as the proof of the corresponding estimates for and S in [JNRZl], noting that 
the substitution of factor (0(^) — (f)(0) — ^d^(f){0)) in R for factor (^(C) — (t){0)) in S introduces an 
additional factor of (1 + in the decay rate for R. □ 

We require also the following key new estimates, proved by similar techniques. The first statement 
reflects the fact that modulations approximately travel at reference group speed a. The second one 
quantifies the fact that the (parabolic) second-order linearized modulation equation is given by 
(3.13). 

Lemma 3.3. Under assumptions (H1)-(H2) and (D1)-(D2), for allt>0, 2 < p < oo, 
(3.10) 11(5* - a5,)(.^'(t)5)|L.(M) < min |5J + ft'~'J_'!~'^^^^ 



(3.11) ||(9t-a9,)(sf(t)(W))llL.(R) < (l + i)-^('-'/^)-^||9./io||Li(M)- 
Proof. Differentiating (3.2), we have 

(3.12) {dt - ad,){sP{t)g){x,t) = T (A(0 - aiOe'iMOe^^^^'{m,m)L^i[o,i]){m 

J —IT 

with A(^) — ia^ = 0(|^p), whence the result follows, again by (2.4). □ 

Proposition 3.4. Assuming (H1)-(H2) and (D1)-(D2), let a{t) he the solution operator of the 
convected heat equation 

(3.13) ut = aux + duxx, 

where a, d are as in (2.5), and let f be a periodic function on [0, 1], / G L^(]R). Then, for all t > Q, 
2<p<oo, and I G N, 

(3.14) \\disP{t){hof) - imJ) <^md%)\\LviR) < {l + t)--^^'-'/Ph--2\\dxho\\mR)nLHRy 
Proof. Expressing 

(3.15) (0(0), /)i2([o,i]) Ht)d%){x) = [ (iO^e^^-e(-«-'^«')*(^(0), /)i2([o,i]) hoiOd^ 
and (recalling (2.2)) 

J-n 



and subtracting, we obtain, for I G N, 
[(^(0), f)a{t){dlho) - diisPmfho))] (x) = / e^^-e(-^-'^«')*(zO'-^^(0), /) djToiOd^ 

J — TT 
J —IT 
J —IT 

each term of which may readily be estimated similarly as in the estimates on s'^ in [JNRZl], 
using (2.4), the fact (by (D2)) that KA(0 < > 0, the Hausdorff- Young bound ||dA)||L- < 

lia^/iolUi, and the estimates A(0-(i<-d^') = O(|eP),<^(0-^(0) = 0(|ei), and (l-a(O) < |?| to 
yield the result. In the estimate of the last term on the righthand side, we are using /e^*^^') = 

{^{Of)'ji where 4>{0fj denotes the jth coefficient in the Fourier expansion of periodic function 
0(0/ 5 to estimate 



<\\d.ho\\LKY.MOf)j\' 



< 



V j¥o 

ll^(0/llL2([0,l])l|5a:/io||Li(K)- 



□ 



4. Nonlinear stability estimates 



4.1. Perturbation equations and nonlinear decomposition. Refining [JNRZl], for u sat- 
isfying ki,ut = k^Uxx + f{u) + ki,cux and ip{x,t) determined as in Proposition 1.1, we set first 
v{x, t) = u{x — ip{x, t),t) — u{x) and then 

(4.1) z{x, t) = u{x — ip{x, t), t) — u{x) — dku{x)ki,dxip{x, t) = v{x, t) — dku{x)k-),dxip{x, t). 

Also we set do := uoi- — ho{-)) — u and introduce % a smooth cutoff such that x(t) is zero for f < 1/2 
and one for t > 1. 

Lemma 4.1. The nonlinear residual z defined in (4.1) satisfies 

z{t) = R{t){do + hou') + [ R{t-s)M{s)ds 

Jo 

+ (1 - x(i)) (BP{t){do + hou') - hou' - dkU kAho + / RP{t - s)Ar{s)ds) , 

(4.2) ^ 

tp{t) = sP{t){do + hou') + [ sP{t - s)Mis)ds 

Jo 

- (1 - x(t)) (sP{t){do + hou') -ho + J^ sP{t - s)Af{s)ds^ , 



where 

(4.3) KM =Q + n^ + {Kdt + kldl)S + T 
with 

(4.4) Q := f{v + u)- f{u) - df{u)v, 

(4.5) 7^ := -Kvi)t - klvtpxx + kl{ux + v^) ^ , 

(4.6) S := v^j^, 
and 

(4.7) r:=-{f{v + u)-f{u))^,. 

Proof. Recall from [JNRZl] that we have chosen ip so that i/j and v satisfy both (4.2) (ii) and 

v{t) = S{t){do + hou') + [ S{t- s)M{s)ds 

(4.8) 

+ 







(1 - xit)) (u'sP{t){do + hou') - hou' + u'sP{t - s)Mis)ds) , 

where S{t) = S{t) - u'sP{t). Using now z{t) = v{t) - dkuk^^x{t), R{t) = S{t) - dkukAs^{t) and 
RP(t) = {u' + dkuk^dx)sP{t), equation (4.2)(i) follows from (4.2)(ii) and (4.8). □ 

4.2. Refined stability estimate. With these improvements, we may refine Proposition 1.1. 

Proposition 4.2. Under the assumptions of Proposition 1.1, for t > 0, 2 < p < oo, we have the 
sharpened estimate 

(4-9) \\z{t)\\LrQ&)<Eoln{2 + t) {l + tyl 

Proof. Prom the proof of Proposition 1.1 in [JNRZl] ,^ we have the crude bounds, for 2 < p < oo, 

II^(*)IIlp(m), m,i^xm\\w'<+^,nm ^ CEo{i + t)-"^^'-'/p^ 

(4.10) J 

\\vmHK(^^)<CEo{l + t)-^, 

whence 

(4.11) \\MmH^M)<cEoii+tri 

Applying the bounds of Propositions 3.1 and 3.2 to system (4.2), we obtain for any 2 < p < oo 

Mt)\\Lpm < CEoil + i)-^(i-Vp)-| +CEo {l + t- s)-^(V2-i/p)-1(i + ^yUs 

(4.12) Jo 

< CEoln{2 + t) {1 + tyi, 
verifying the result. □ 



"'Though for simphcity only bounds were reported in the Proposition, the stated higher-derivative bounds were 
obtained in the course of the proof. 



5. The Whitham equation 
We complete our investigation by connecting to the Whitham modulation equation. 
Lemma 5.1. Assuming (H1)-(H2), (D1)-(D2), KM{t) = P kltpxitf + r{t), where 

(5.1) = f {u){dku, dku) + Kdku" - -^df{u)dku + u" - adkv! 
is "periodic and \\r{t)\\ii(^-^ < Eq{1 + t)~^. 

Proof. Immediate from Lemma 4.1, (4.10), (4.9), and ||V't(t) - ailJx{t)\\L^ = 0((1 + t)"^/^), a 
consequence of Lemma 3.3. □ 

Lemma 5.2. Setting e := (^(0),/^), we have e = ^^"(A;*). 

Proof. We may rewrite 

fP = -LdkU + 2Kdku" - J{K)dku' + u" + ]^d^f{u){dku, dku) 

thus 

e = ((^(0), 2Kdku" - u'{K)dku' + u" + ]^d'' f{u){dkU, dku)) = . 

The last equality comes differentiating twice the profile equation with respect to k (see [DSSS]). □ 

Proof of Theorem 1.3. With estimates on k^^dxi^it) from Proposition 1.1, the estimate on u follows 
from 

(5.2) ||u(--V'(-,t),t)-^/'*(-)-5fc«''*(-)MxV^(-,t)||LP(R) <i?oln(2 + t) {l + tri 

established in Proposition 4.2 for 2 < p < oo. 
Using Duhamel's principle we may write 

(5.3) k{t) = a{t)kAho + [ <j{t- s)d^{ek\s)))ds , 

JO 

where a is the constant-coefficient solution operator defined in Proposition 3.4. On the other hand, 
by Propositions 3.1, 3.2, and 3.4, 

ft 



(5.4) 



k^'^x{t) = k^x{t) (^dxsP{t){do + hou') + ^ dxsP{t - s)M{s)ds^ + (1 - x{t))kAho, 

= a{t)kAh{i + / (j{t- s)dx{e{k^i)xf{s))ds + f{t), 
Jo 

where 

\\m\\Lp(R)<EoH2 + t) (l + 0-5(^-^/f)-^, 2<p<oo. 
(Here, we have used (3.14) and 

\\Mi^lis))\\L^R)nL^R) ^ \\tpxis)\\L^R)nL'>o(R)\\tpxxis)\\L^R) ^ (1 + s)"^ 

to bound 

r dxsP{t - s)fPkl^x{tfds- f a{t - s)dx{e{k^^xf{s))ds\\L,(^) 
Jo 

Jo 

pi 

< / {i+t-s)^'^^-^/p^{t-s)-ki+sr'ds.) 

Jo 



Thus, subtracting (5.3) from (5.4), and defining S := ki.ipx — k, we have 



(5.5) 6it) 
Defining, for some i] > 0, 



ait - s)dx{e6{s)ik{s) + k^i;^is)))ds + r{t). 



u{t):= sup sup ||5(s)|U.(K)(l + s)5(i-Vf)+i-'7, 

p6[2,oo] 0<s<t 



we thus obtain by the standard heat bounds ||<7(t)5a;/||x,p( 
1 < 9 < P < oo, that, for 2 < p < oo, 

t/2 



< cr^(v?-i/P)-| 



\ds 



when 



WmWLviM.) < I (t- s)--2^'-'M-'2\\e6is)ikis) + k^MsMLH 
Jo 

+ f\t- s)-5(V2-i/f)-^||e5(5)(A;(5) + k,MsmL^iR)ds + wrnhn 

Jt/2 
rt/2 

Jo 

+ / {t-s)-^^^/^-^/P^-^i^{t)Eo{l + s)-^+''ds 

Jt/2 



+ Eo ln(2 + t) (1 + t)-5(i-i/p)-| 

<Eo(K^) + l)(l + ^)-^('-'/^^-^+^ 



giving 



Uit) < Cr,Eo (1 + Uit)) , 

whence, if Eq < 1/(2C^), we have (noting that u by standard theory remains bounded) that 

iy{t) < 2Cr,Eo. This gives ||A;(t) - A;,5,V'(t)||Lp(M) = l|5(i)llLp(R) < 2C^(1 + t)-5(i-Vp)-'?, completing 
the result for k — /c^Vx- A similar computation yields the result for /i — ^. □ 

Acknowledgement. Thanks to Bjorn Sandstede for pointing out the results of [SSSU]. 
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Appendix A. Asymptotic equivalence of scalar equations 

Lemma 1.2 follows by a simplified version of the arguments used in [Ka, LZ] to prove correspond- 
ing but weaker versions in the system case. We include a proof here, both for completeness, and to 
motivate the more complicated comparison arguments appearing in the main body of the paper. 

Proof of Lemma 1.2. By the general, system, results of [Ka],^ provided Eq := ||fco||Lin//3(R) is 
sufficiently small, we have for := k — fc* and 1 < p < oo, the "heat- type" bounds 

^ ' ~ _3 

llfcccWllffl(R), ||A;a;(t)||ifl(R) < £^o(l+i) ^■ 

Setting (5 := K — K, we have, subtracting and rearranging, 

(5t - a4 - db^^ = d^T, T = 0{i\k\ + \~k\)5) + 0(fc3) + oCkk), 
with S\t=o = 0, where a = dg(A;*), d = d{ki,) are constant. By Duhamel's formula, 

S{t) = [ a{t- s)d^J^{s)ds, 
Jo 

where a is the solution operator of the convected heat equation ut — aux — duxx = 0. Applying the 
standard heat bounds ||<7(i)5^/||iP(iR) < t 2*-? p> 2 1 < q < p < 00, together with 



See [Ka, Remark 4.2] improving the result of [Ka, Theorem 4.2] in the strictly parabolic case. Bounds (A.l) are 
proved by estimates similar to those of this paper and [JNRZl]; see for example (A. 2) just below. 
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(A.2) t^l"^ 



< 



3 

1 < 9 < 2, we find, defining vit) := supQ<5<( ||5(s)||i2(]R)(l + s)i~^, that, for all 1 < p < oo, 

rt/2 

Jo 

rt/2 

Jo 

Jt/2 

whence < C^^Jq (£^o + '^(O)- This implies that u{t) < 2CnE^ for Eq < 1/(2C^), giving 

< 2C,i?2(i 1 < p < oo. 

Finally, let 4i{x,t) = (^(x/\/F+T) define a self-similar solution of Burgers equation (1.16) such 
that J (p = J ko. Defining 6 := k — cj), we have by the and first moment assumptions that, 
for 1 < p < DO, 

Expressing by Duhamel's formula 

5[t) = a{t)6{0) + [ a{t- s)d^f{s)ds, 
Jo 

where T = 0{{\k\ + It/'DfJ), and, for any ij' > 0, estimating as before, we obtain by a contraction 
argument similar to the above that 

WmiLPiR) < E, (l+i)-|(l-l/p)-|+V, 1<P<00, 

provided Ei is small enough (depending on r]'). Indeed one may prove for 

i>it):= sup sup ||5(s)|U.(M)(l + s)i(i-Vp)+^V 

pG[l,oo] 0<s<t 

that £>(t) < C'^,Ei{l + i>{t)) and conclude that i>{t) < 2C[^,Ei if C[^,Ei < 1/2. □ 

Remark A.l. The order 0{dx{k6)) ~ 0{dx{kkx)) ~ 0{dx{k^)) of neglected terms in the argument 
above is consistent with the order neglected throughout the paper. Indeed, as emphasized by 
Remark 1.4, the true local wavenumber is A;*^'^ (where ^(^t) is the inverse of y i— >■ X{y,t) := 

y — ilj{y,t)) which differs from A;*(l + ipx) by a 0(|Ki'x|P) term since 

^x{x, t)-[i + Mx, t)] = M^ix, t),t)- M^{x, t) - ^(*(x, t),t),t)^x{x, t){i - M^{x, t), t)). 

With a slight bit of additional effort, one can check that k^^^x satisfies the Whitham equation or 
its quadratic approximant only up to a truncation of the order 9^(V'x) ~ dx{kkx)- Once integrated, 
this means that the formal Whitham equation (1.13) derived for fc^(l + V'a;) can be expected to hold 
for ki,^x only to the (asymptotically equivalent) quadratic order considered here. In particular, if 
wc were to derive higher order expansions for the equation satisfied by ^x^ as we could in principle 
do, these would not in general agree with the Whitham equation (1.13). In any case, while equation 
(1.16), most easily obtained by switching to a comoving reference frame x ^ a; — at, is designed to 
match, about a given wave, exactly the order of description we attain here, equation (1.13) is just 
a convenient way to piece together the two first orders of a nonlinear WKB expansion; thus (1.13) 
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could in principle be meaningful even when dealing with a modulated background wave (instead 
of a true wave) but should not improve in accuracy its quadratic approximant about a given wave. 
For further discussion of this and related issues, see [DSSS, NRl, NR2]. 
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